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DTIC  QUALTIT 


EXTENSION  OF  FRESNEL  NEAR-FIELD  REGION 
BY  SINGLE-CYCLE  SINUSOIDAL  PULSES 

1.  INTRODUCTION 

The  theory  of  diffraction  of  a  monochromatic  (or  a  CW)  wave  [1,2]  by  a  circular  disk  (or  aper¬ 
ture)  of  radius  a—DH,  shows  that  the  diffracted  energy  density  in  the  Fresnel  near-field  region 
decreases  slower  than  1/z^  on  the  axis  of  the  disk.  Beyond  this  region  the  well-known  behavior  is 
l/z^  in  the  Fraunhofer  region  (far  field  region).  Though  there  is  no  well-denned  demarcation 
between  the  Fresnel  and  the  Fraunhofer  region,  the  optical  [1,2]  and  the  acoustical  [3]  literatures  take 
this  to  be  /X,  where  X  is  the  wavelength  of  the  wave  concerned.  On  the  otherhand,  it  is  a  common 
practice  in  the  antenna  literature  [4]  to  assume  this  boundary  to  lie  between  D^/X  (=  4o^/X)  and 
2£)^/X  (=  8d^/X).  An  approximate  theory  [1]  can  be  used  to  explain  that  as  the  point  of  observation 
z  from  the  disk  (or  aperture)  is  made  longer  than  n^/X,  there  is  less  than  one  Fresnel  zone  in  the  disk 
and  the  intensity  of  the  wave  energy  decays  as  1/z^  (as  in  the  Fraunhofer  region  or  far  field).  On  the 
otherhand,  as  z  is  made  to  shrink  towards  the  aperture  along  its  axis,  so  that  more  than  one  Fresnel 
zone  in  the  aperture  contributes  to  the  intensity  at  z,  the  energy  density  oscillates  between  zero  and  a 
maximum  value.  For  extremely  small  z,  this  approximate  theory  is  not  applicable.  The  first  max¬ 
imum  intensity  appears  at  z  =  a^/X,  beyond  which  the  wave  intensity  starts  to  decrease  like  1/z^. 
This  observation  then  indicates  that  the  Fresnel  near-fleld  region,  where  the  energy  density  drops 
slower  than  1/z^,  increases  as  the  frequency  of  the  incident  field  increases.  The  decrease  of  intensity 
slower  than  1/z^  implies  also  that  the  diffracted  or  radiated  field  is  collimated. 

In  many  applications  which  require  transmission  of  high  field  intensity  (such  as  an  extended 
radar  system),  to  a  relatively  long  distance  (a  few  hundred  to  a  few  thousand  kilometers),  in  principle 
this  distance  can  be  made  to  fall  within  the  Fresnel  near-fleld  region  by  simply  increasing  the  fre¬ 
quency  of  the  field  as  well  as  the  aperture  dimension  of  the  radiator.  However,  just  increasing  the 
frequency  of  the  monochromatic  source  of  the  radiator  may  not  be  an  efficient  method,  as  indicated  in 
the  following  recent  works  [5,6,7, 8]. 
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Wu  [5]  who  introduced  the  concept  of  "eiectromagnetic  missile*,  defining  the  electromagnetic 
energy  which  decays  slower  than  1/r^,  showed  that  some  pulsed  sources  backed  by  a  reflector  can 
transmit  electromagnetic  energy  to  a  great  distance,  where  the  field  energy  falls  slower  than  the 
inverse  of  the  square  of  the  distance.  Wu  and  Shen  [6]  also  preformed  an  experiment  in  a  laboratory 
for  demonstrating  the  validity  of  the  theory.  Besides  them,  Ziolkowski  and  his  associates  [7,8]  also 
investigated  the  behavior  of  a  special  type  of  highly  collimated  wave,  named  localized  wave 
transmission*  which  also  required  an  excitation  of  an  aperture  by  suitably  created  pulses.  They  [8] 
also  conducted  an  acoustic  experiment  for  verilying  their  concept  and  theory.  The  aims  of  these  two 
groups  ([5,6],  [7,8])  appear  to  be  the  same,  though  their  respective  theoretical  approaches  and  con¬ 
cepts  are  different.  For  instance,  in  the  study  of  'electromagnetic  missiles,*  a  conventional  method 
of  solution  of  wave  equation  in  the  Fresnel  region  is  needed,  which  makes  it  easier  to  understand  the 
process  of  radiation.  However,  a  special  solution  (appears  to  be  unconventional)  together  with  the 
choice  of  a  special  type  of  pulsed  sources  are  required  in  Ziolkowski’ s  is  work.  In  spite  of  the 
apparent  conceptual  and  methodological  differences,  the  studies  of  both  the  'electromagnetic  missiles* 
and  the  *localized  wave  transmission*  show  that  directed  short  pulses  encompassing  a  wide  band  of 
high  frequencies  are  required.  The  results  of  these  studies  provide  a  clue  to  the  question  of  finding  a 
method  of  extending  the  Fresnel  near-field  region  more  efficiently. 

In  references  [5,6]  non-sinusoidal  pulses  with  suitable  forms  were  used  so  that  the  results  could 
be  expressed  and  interpreted  conveniently.  For  non-sinusoidal  pulses  there  is  no  unique  and  apparent 
equivalent  wavelength,  corresponding  to  which  a  Fresnel  region  can  be  identified.  As  a  result,  there 
is  no  convenient  and  well  accepted  way  to  determine  the  enhancement  of  the  near-field  region,  when 
the  source  of  radiation  is  a  non-sinusoidal  pulse.  On  the  otherhand,  for  a  sinusoidal  pulse,  such  as  a 
single-cycle  sinwot  or  cosoiot,  the  extension  of  the  Fresnel  region  can  be  ascertained  in  a  convenient 
manner  by  comparing  the  respective  results  with  those  of  a  CW  sinojof  source. 
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In  view  of  the  above  observations,  the  present  investigation,  resembling  the  study  of  "elec¬ 
tromagnetic  missiles, '  studies  the  feasibility  of  extending  the  Fresnel  near-field  region  of  a  conducting 
circular  disk  antenna  [Fig.  1]  excited  by  sinusoidal  pulses,  single-cycle  sincoot  and  single-cycle  cosojqI, 
which  are  more  familiar.  It  may  not  be  too  difficult  to  generate  a  short  pulse  containing  a  few  cycles 
of  sincdot.  The  consideration  of  a  single-cycle  cosuQt  appears  to  be  more  academic  than  realistic, 
nevertheless  the  results  shed  some  light,  pointing  out  clearly  that  the  shorter  the  rise  time  of  a  pulse, 
the  slower  is  the  decay  of  the  radiated  energy  density.  The  efficacy  of  a  single-cycle  sinuot  pulse 
over  the  corresponding  CW  sincjot  as  exciting  sources  in  extending  the  Fresnel  near-field  region  (and 
hence  collimating  the  radiated  energy)  lies  in  the  cohesive  contributions  of  all  the  higher  frequencies 
equal  to  and  above  ci)o  (i.e.  all  oi  in  o)  >  oio)  contained  in  the  pulse. 

In  section  2  presentation  and  discussion  of  numerical  results  are  made,  followed  by  conclusions 
in  section  3.  Development  of  the  theory  of  radiation  in  the  Fresnel  near-field  region  of  a  circular 
conducting  disk  excited  by  uniform  sinusoidal  current  pulses  and  the  associate  analyses  can  be  found 
in  Appendixes  A  and  B. 

2.  DISCUSSIONS  OF  NUMERICAL  RESULTS 

Since  the  behavior  of  the  energy  density  as  a  function  of  the  coordinate  z  along  the  axis  of  the 
conducting  disk  or  the  dimensionless  distance  a(=z/(D^/Xo))  will  demonstrate  whether  these 
sinusoidal  pulses  do  indeed  enhance  the  Fresnel  near-field  region,  we  shall  first  present  the  energy 
densities  for  different  parameters.  The  wavelength  X©  defined  by  Xq  =  Ire /u>o,  where  c  is  the 
velocity  of  light  in  free  space.  In  Fig.  3a  three  normalized  energy  density  expressions  associated  with 
a  CW  sinuot  current  (Eq.  24,  Appendix  A),  a  single-cycle  sinwot  pulse  (Eq.  34a,  Appendix  A)  and  a 
single-cycle  coswot  (Eq.  34b,  Appendix  A)  are  presented  as  functions  of  a  along  the  axis  of  the  disk. 
The  energy  per  cycle  of  each  of  these  excitations  is  the  same. 
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The  behavior  of  the  radiated  normalized  energy  density,  Pcw/Pq  =  sin^ 
to  the  CW  sincdot  current  source  is  well-known  [3,5,6].  It  is  highly  oscillatory  for  a  :£  1/8,  reaches 
the  maximum  value  1  at  a  =  1/4,  and  then  begins  to  drop  monotonically.  This  oscillatory  behavior 
of  Pew  for  a<  1/8  is  due  to  the  fact  that  more  than  one  Fresnel  zone  in  the  disk  (source  region)  con¬ 
tributes  to  the  illummation  and  causes  interference.  Here  we  make  a  distinction  between  the  Fresnel 
zone  and  the  Fresnel  region.  The  Fresnel  zone  is  located  in  the  source  region,  whereas  the  Fresnel 
region  lies  in  front  of  the  radiator,  where  near-field  radiation  takes  place. 


In  contrast  the  normalized  radiated  energy  density  Ps/Pq,  associated  with  the  single-cycle  of 
simoQt,  has  a  constant  value  unity  without  showing  any  oscillatory  behavior  in  the  region  a  <  1/8. 
For  a  >  1/8,  it  increases  to  a  maximum  value  of  1.5  at  a  =  1/4,  from  there  it  decreases  monotoni¬ 
cally  with  the  increase  of  a.  The  oscillatory  behavior  of  P^w  for  a  <  1/8  is  also  absent  in  Pc/Pq. 
the  normalized  energy  density  corresponding  to  the  single-cycle  costoot  excitation.  This  observation 
suggests  that  if  there  are  nulls  in  the  radiation  pattern  of  a  radiator  excited  by  a  monochromatic  sig¬ 
nal,  those  nulls  can  be  eliminated  by  choosing  a  short  pulse,  the  spectrum  of  which  conuins  that 
monochromatic  source  frequency  as  a  source  of  excitation. 

However,  the  behavior  of  Pq/Pq  shows  a  peculiarity,  since  it  begins  to  drop  from  unity  to  a 
minimum  value  (about  0.78)  near  a  =  0.162,  from  where  it  starts  to  increase  to  a  maximum  value 
(about  1.58)  at  a  =  0.3.  Before  attaining  this  maximum  value  Pc/Pq  passes  through  the  point  of 
maximum  of  P,/Po  at  a  =  1/4.  From  a  =  0.3,  Pc/Po  begins  to  decrease  monotonically  with  q. 
although  at  a  slower  rate  than  either  PJPq  or  Pew^Po^  exhibiting  the  fact  that  P^  has  a  larger  value 
than  that  of  P,  or  P^w  .  Between  o  =  1/8  and  1/4,  P^  is  smaller  than  Pj.  However,  both  P,  and  P, 
are  larger  everywhere  than  Pew,  if  we  agree  that  the  average  value  of  Po,./Po  is  only  1/2  for 
a  <  1/8,  since  both  P,/Po  and  Pf/Po  are  unity.  Although  some  graphs  are  shown  for  a  =  0  (or 
z  =  0),  they  are  not  valid.  However,  a  can  be  very  small,  depending  on  the  value  of  Xq  and  the  con¬ 
dition  a  /z  «  1 . 
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The  larger  values  of  and  Pc  are  due  to  the  simultanraus  and  collective  contributions  of  all  the 
high  fi«quencies  exceeding  coq  (w  ^  uo)  contained  in  the  single-cycle  pulses.  It  is  not,  however,  clear 
why  Pc/Pq  is  smaller  than  P,/Po  in  the  region  1/8  <  a  <  1/4.  It  may  be  that  for  a  certain  large 
the  amplitudes  of  the  frequency  spectrum  between  uq  ^  u)  S  (i)„,  of  the  single-cycle  cosuqI  are 
smaller  than  those  for  the  single-cycle  sin  ojQt,  although  the  single-cycle  cos  uqI  contains  more  high 
frequencies  with  large  anq)litudes  in  general.  Note  that  the  spectra  of  the  single-cycle  sin  oiqI  and  the 
single-cycle  cos  UQt  behave  like  l/u^  and  l/oj,  respectively,  as  w  —  oo. 

Figure  3a  shows  the  behavior  of  P^,  P,  and  P^  for  small  values  of  a(<  1).  On  the  other  hand. 
Fig.  3b  displays  how  P^h,,  Pj  and  P^  decrease  for  larger  values  of  a,  which  is  helpful  in  determining 
whether  P,  and  P^  have  indeed  shown  extension  of  the  near-field  region  when  compared  with  Pew- 
For  this  purpose  let  us  choose  first  a  value  of  a,  which  we  take  to  be  2  (i.e.  z  =  2D^/\q)  for  con¬ 
venience.  At  this  distance  along  the  axis  of  the  disk  Pew/^o  ^  value  about  0.03806,  which  is 
equivalent  to  -14.2  dB.  We  now  find  from  Fig.  3b  that  P^/Po  and  Pc/Pq  drop  to  this  value  at 
a  —  11  and  30  respectively.  This  finding  can  be  interpreted  as  follows:  when  compared  with  a  CW 
sin  ojot  excitation,  the  single-cycle  sin  uo^-and  the  single-cycle  cos  upt  enhance  the  Fresnel  near-field 
region  by  37.5%  and  275%,  respectively.  Although  an  ideal  single-cycle  cosupt  is  not  realizable 
(similar  to  the  case  of  an  ideal  rectangular  pulse,  which  is  used  in  many  examples  in  the  literature  for 
the  purpose  of  illustrations)  the  associated  results  shed  light  in  showing  clearly  that  the  higher  the  rate 
of  rise  of  an  exciting  pulse,  the  slower  is  the  decay  of  the  radiated  energy. 

Let  us  now  mm  our  attention  for  a  moment  to  an  interesting  aspect  of  the  initial  time  variation 
of  a  pulse  and  the  corresponding  behavior  of  its  spectrum  at  high  frequencies.  Reference  [6]  shows 
that  if  the  current  pulse  behaves  like  for  small  values  of  t,  then  its  spectrom  P(«)  at  high  fre¬ 
quencies  decreases  like  a)~^‘  Consequently,  the  corresponding  radiated  energy  density  P(z) 

along  the  axis  of  the  disk  drops  as  z  with  increasing  z-  In  order  for  a  current  pulse  to  produce  a 
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"missile*  effect,  they  required  0  <  e  <  1,  so  that  the  energy  density  decreases  slower  than  1/z^. 
For  a  single-cycle  sinci>ot  we  can  not  apply  these  results  in  a  strict  sense.  For  instance,  for  small  t, 
siouQt  increases  as  t.  If  we  use  their  results,  we  get  e  =  3/2,  showing  F,{ui)  -  1/u^  and 
Ps/Pq  —  l/z  ,  without  displaying  any  "missile*  effect.  Although  the  behavior  of  Fj(w)  is  correct  as 
<i)  —  00,  our  result  [Eq.  34a,  Appendix  A]  shows  that  Ps(z)/Po  behaves  like  1/z^  for  very  large  z. 
In  spite  of  this  behavior  of  P^iz),  it  displays  [Figs.  3a  and  3b]  effective  enhancement  of  the  near¬ 
field,  i.e.  showing  "missile"  effect.  Using  a  similar  argument  for  the  single-cycle  costoot,  one  finds 

A 

that  for  small  /,  costijo?  =  1  indicating  e  =  1/2,  Ff(«)  ~  \/o>  and  Pc  ~  l/z  for  large  z.  Surpris¬ 
ingly,  such  behaviors  agree  with  our  results,  although  the  single-cycle  cosuiQt  is  unrealizable. 

The  above  discussions  pertain  to  the  axial  behavior  of  the  energy  densities,  which  are  not  too 
difficult  to  compute.  On  the  other  hand,  the  off-axis  (i.e.  s  =  p/a  >  0)  behavior  of  either  P,  or  Pc 
is  very  difficult  to  compute,  since  the  expressions  (32a)  and  (32b)  of  Appendix  A  converge  slowly. 
In  f  ~  p/a,  p  is  the  perpendicular  distance  from  the  z-axis.  However,  for  s  =  1  and  j  » 1,  it  is 
relatively  easy  to  calculate  P,  and  Pc  in  closed  form  [Eqs.  (39a)  to  40b,  Appendix  A].  For  s  » 1 
the  associated  infinite  series  can  be  adequately  approximated  by  retaining  only  the  first  term.  In  Figs. 
(4a)  and  (4b)  we  show  these  results.  Figure  4a  displays  the  behavior  of  the  normalized  energy  den¬ 
sity  Ps/Pq  for  5=  0,  0.125,  0.25  and  1.  Similarly,  Fig.  4b  presents  Pc/Pq  for  the  same  values  of 
5.  The  cases,  s  =  0.125  and  0.25  for  Ps/Pq  are  computed  by  using  Eqs.  (36a)  and  (36b)  of  Appen¬ 
dix  A,  respectively.  For  s  »1,  both  normalized  energy  densities  become  so  small,  that  they  are 
excluded  from  these  displays.  In  general  these  results  show  that  both  P^  and  Pc  are  concentrated  in 
the  vicinity  of  the  disk  axis  and  drop  rapidly  with  increasing  values  of  s.  Such  behaviors  of  Pj  and 
Pc,  together  with  an  examination  of  the  expressions  (24)  to  (27)  of  Appendix  A,  indicate  that  the  radi¬ 
ated  energy  density  associated  with  a  short-pulse  excitation  is  highly  focused  towards  the  axis  of  sym¬ 
metry  of  the  radiator  and  extends  to  a  larger  distance  along  this  axis  when  compared  to  the 
corresponding  results  for  the  CW  excitation. 
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In  Figs.  Sa  through  Sd  the  behavior  of  the  normalized  radiated  electric  field,  E„,  associated 
with  the  single-cycle  sinuot  are  presented.  The  corresponding  behaviors  of  E„  (oo*  shown)  associated 
with  the  single-cycle  coscd^r  can  be  obtained  by  taking  the  first  derivative  of  £„  with  respect  to 
{t*  =  t  —  z/c).  Figures  Sa  through  Sd  show  as  a  function  of  the  normalized  time  t  =  t  n, 
where  T  =  Iv/iOQ  is  the  period  of  one  cycle,  for  r  =0  and  1  at  a  given  distance  a,  parallel  to  the 
z-axis.  The  thick  dotted  curve  is  common  to  both  ^  =  0  and  1.  It  follows  from  Eq.  28a  of  Appendix 
A  that  for  s  =  0,  the  two  signals,  one  radiated  from  the  center  of  the  disk  [the  first  term  of  (28a)] 
and  those  radiated  from  the  circumference  of  the  disk  (second  term),  overlap  for  a  <  1/8.  For 
s  =  1  [Eq.  30a,  Appendix  A],  the  radiated  fields  from  these  two  regions  overlap  for  a  >  1,  other¬ 
wise  they  are  separated.  Thus,  Figs.  Sa  through  Sc  show  that  these  two  radiated  signals  overlap  for 
s  =  0,  and  are  well  separated  for  r  =  1.  However,  in  Fig.  Sb  one  finds  that  for  r  =  1,  the  two  sig¬ 
nals  begin  to  combine.  The  second  signal  (s  =  1),  which  comes  from  the  center  as  well  as  from 
other  circumferential  points  that  are  farther  away,  is  much  smaller  than  the  first  signal,  which  comes 
only  from  the  nearest  point  on  the  circumference.  The  interference  of  these  rays  among  themselves 
reduces  the  amplitude  of  the  second  signal.  Figure  Sd  shows  that  for  a  large  s(=5),  the  off-axis  field 
strength  is  very  small,  which  agrees  with  what  is  intuitively  expected. 

3.  CONCLUSIONS 

It  has  been  shown  that  in  order  to  extend  the  Fresnel’s  near  field  region,  which  process  being 
the  same  as  collimating  (or  focusing)  the  radiated  energy,  infront  of  a  conducting  circular  disk 
antenna,  a  single-cycle  sinoiot  current  pulse  is  more  efficient  than  the  corresponding  CW  sinuot 
current  source.  Such  an  ability  of  a  single-cycle  pulse  lies  in  the  cohesive  contributions  of  the  high 
frequencies  above  (Oq  (o’  2  ciio)  contained  in  the  pulse.  Assuming  that  the  near  field  region  extends 
to  unit  distance  (along  the  disk  axis)  due  to  the  excitation  by  a  CW  sinoior  current,  it  is  found  that  the 
corresponding  distances  are  1.375  and  3.75  when  the  exciting  current  sources  are  a  single-cycle 
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siocdot  and  a  single-cycle  coswoi,  respectively.  The  use  the  single-cycle  cosu)Qt  pulse  may  be  regarded 
as  OKHre  academic  than  realistic,  nevertheless  the  corresponding  results  show  clearly  that  the  shorter 
the  rise  time  of  a  pulse  consisting  of  large  amplitudes  of  high  frequency  components,  the  slower  is  the 
rate  of  fall  of  the  radiating  energy,  thereby  proving  itself  more  efficient  in  focusing  or  collimating  the 
field  energy  in  the  vicinity  of  the  axis  of  symmetry  of  the  radiator.  If  there  is  a  way  to  generate  a 
short  pulse,  the  hequency  spectrum  of  which  falls  like  \/b/  as  u  increases  indefinitely,  with  0  lying 
at  least  in  the  range  1  s  /3  <  2,  it  should  be  welcomed  as  a  potential  candidate  for  the  applications 
where  the  above  mentioned  behavior  of  the  radiated  energy  in  the  Fresnel  region  is  desirable. 

The  results  of  this  study  indicate  also  that  if  there  are  nulls  in  the  radiation  pattern  of  a  radiating 
system  excited  by  a  monochromatic  signal,  those  nulls  can  be  made  to  disappear  if,  instead,  a  short 
pulse,  the  spectrum  of  which  contains  that  frequency  of  the  monochromatic  source,  is  used. 
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APPENDIX  A 


I.  GENERAL  FORMULATION 

For  the  transmission  of  a  directed  and  collimated  electromagnetic  energy  one  needs  a  focusing 
device  in  addition  to  an  antenna.  However,  for  the  purpose  of  demonstration  and  simplicity  we  shall 
take  a  perfectly  conducting  circular  disk  of  radius  a  as  an  antenna,  which  also  serves  as  a  focusing 
mechanism.  This  circular  disk  antenna  is  excited  by  sinusoidal  single-cycle  uniform  current  pulses, 
sinojo/  and  cos  uqI,  respectively.  The  primary  reason  of  our  choice  of  these  sinusoidal  pulses  is  to 
show,  by  comparing  with  the  result  corresponding  to  a  cw  sinusoidal  excitation,  whether  these  pulses 
do  indeed  extend  the  Fresnel’s  near  field  region. 

At  first,  the  problem  will  be  formulated  in  the  frequency  domain,  and  then  the  corresponding 
time  domain  result  will  be  obtained  by  an  application  of  the  inverse  Fourier  transform.  It  should  be 
noted  that  the  assumption  of  exciting  the  circular  disk  by  a  uniform  current  imposes  an  unnecessary 
restriction,  which  implies  filtering  out  very  .high  fr«]uencies.  Nevertheless,  it  is  adequate  [S,6]  for 
the  purpose  at  hand.  Consider  a  thin  perfectly  conducting  disk  of  radius  a,  lying  in  the  x-y  plane  in 
otherwise  free  space  [Fig.  1].  The  center  of  the  disk  coincides  with  the  origin  of  the  coordinate  sys¬ 
tem.  This  disk  anteiuia  is  excited  by  a  uniform  x-directed  surface  current  /((r',u),  where  r'  is  on  the 
disk  and  u  =  2rf  is  the  angular  frequency,  which  will  excite  all  the  rectangular  components  of  the 
electromagnetic  field  except  The  non-zero  field  components  can  be  expressed  in  terms  of  the  x- 

A 

directed  vector  potential  Aj(f,u).  Since  we  are  interested  in  the  transmission  or  radiation  of  elec¬ 
tromagnetic  energy  in  the  z  duection  only,  which  is  perpendicular  to  the  circular  disk,  only  the  com¬ 
ponents  Eji  and  Hy  contribute.  The  x-component  of  the  vector  potential  in  the  frequency  domain  is 
gjven  by 

=  0*o/4ir)f  J  Jx  (r\w)  (e^/R)  dS'  (1) 
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The  integration  is  over  the  surface  of  the  disk,  where  R  =  |7  -  r'  | ,  it  =  u/c  -  2ir/X  and  7(x,y,z) 
and  r'{x',y',0)  are  the  observation  and  source  coordinates,  respectively.  In  the  Fresnel  region,  where 
a  /X  »  1  and  z/a  »  1,  the  in  the  exponential  is  approximated  by  [3-6]. 

R  =  1"^  -  r'l  =  z  -t-  (jf  -  xi^/lz  +  (y  -  y'f/lz  (2a) 

=  z  +  (p^  +  p'^)l2z  -  ipp'/z)  cos  (</>  - 

where 


X  =  p  cos  y  =  p  sin  ^  (2b) 

x'  =  p'  cos  <l>'  and  y'  =  p'  sin  <l>' 

and  p/z  « 1,  a/z  « 1.  In  the  denominator  of  the  integrand  in  (1),  R  is  approximated  by  z.  The 
approximations  indicated  by  (2a)  and  (2b)  with  a /\  >  >1  imply  that  the  Fresnel  region  confines 
itself  close  to  the  axis  of  the  circular  disk.  For  a  uniform  current  let  us  write 

J,  =  /o  ^(«).  (3) 

The  constant  /q  has  a  dimension  ampere/meter.  Then  using  (2a)  and  (2b)  in  (1),  we  have 

The  integral  over  4)'  is  hrJoik  pp'/z),  where  /o(^)  ^  Bessel  function  of  order  zero.  Introducing 

now  the  following  change  of  variables 

f  =  p'/a,  u  -  k  pa/z,  and  7  =  ka^/z,  (5) 

the  expression  (4)  can  be  re-expressed  as 

Ax  (r,a))  =  ^/o  7(w)/(2z)j  exp  Ii^(z  +  p^/2z)]  I{y,u),  (6a) 

where 

I(y,u)  -  Jo  («  f)  exp  {iy  f^/2)  f  </f.  (6b) 
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The  integral  I{y,u)  in  (6b)  cannot  be  evaluated  in  closed  form,  'however,  it  can  be  expressed  in 
terms  of  two  Lonunel’s  functions,  Ui(y,u)  and  U2(y,u)  [9,  10,  11]  in  the  following  manner 

I{y,u)  =  izc/u>)  -  i  (7a) 

where 

Ux  (y,u)  =  y  Jo  (u  n  cos  f7(l  -f^)/2)  r  dt,  (7b) 

=  £  (-ir  (y/u)^^^  J2„^i(u). 

m  =0 

and 

U2(y,u)  =  7  Jo  («f)  sin  (7(1  -  fdf  (7c) 

=  £  (-l)'"(7/«)^^^^2„,  +  2(«)- 

m  =0 

The  representation  of  /(y,u)  given  by  (7a)-(7c),  although  valid  for  all  s  =  p/a  =  11/7,  is 
suited  more  for  large  s.  For  example,  when  s  >  >  1,  only  the  m  =  0  term  is  sufficient  to  obtain  an 
accurate  approximation  of  /(7,«),  i.e.,  I(y,u)  =  (zc/u)Ui(y,u)exp  (iy/2)  =  (zc/w)- 

exp  (iy/2)Ji  (u)/s  +  O  (l/s^\  An  alternative  representation  of  /(y,u),  which  is  suited  for  small 
values  of  5,  can  be  presented  as  follows  [11].  Note  that  the  notations  used  here  are  different  from 
those  in  [11].  Defme  first 


C(7,h)/2  =  Jo(un  cos  (7  f2/2)  fdf. 

(8a) 

S(y,u)/2  =  Jo(un  sin  (7  f2/2)  fdf. 

(8b) 

^'o(^.«)  =  £  (-irs^J2„(u), 

(9a) 

m  “O 

Vx  is,U)  =  £  (-I)”  J2m*x(u). 

m  “0 

(9b) 
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Then  we  have 


Iiy,u)/a^  =  C(y,u)/2  +  iS(y,u)/2  (10) 

C(y,u)  =  (2/y)  [sin  iu^/2y)  +  Vq  (s,u)  sin  (y/2)  -  K,  (jr,«)  cos  (7/2)],  (11a) 

S(y,u)  -  (2/y)  [cos(«^/27)  -  Vq  (s,u)  cos  (y/2)  -  K,  (s,u)  sin  (7/2)].  (11b) 

A  A 

The  field  components  Ejc(r,u>)  and  Hy(f,ui)  can  be  obtained  using  the  following  relations  in  the 
Fresnel  region. 


£,(r,(i))  =  iw  [Aj  +  (l/i^)  Aj,]  =  iu  Ajc(r,u), 


(12a) 


Hy(r,o})  =  (l/no)  Aj,  =  (iu/iiQ)A^  (r,u),  (12b) 

az 

where  jjq  =  V/io/eo  •  Th®  time-dependent  radiated  electric  field  in  the  Fresnel  region  is  then  given 
by  the  following  Fourier  transform. 


Ex(f<0  -  (1/2t)  f  Ex(r,o})  du 

^  — OD 


=  (i  /2ir)  /»o  /o  f  I(y,u)  [«  F(u)/(2z)J  exp  [-i«|r  - 

—  OP 

The  other  quantity  of  interest  is  the  z-directed  energy  density  at  a  point  in  the  Fresnel  region  and 
expressed  as 


Z  +  &- 


/c\]  du 


(13) 


00 

P(r)  =  f  E  (r,t)  X  H(r,t)  •  z  dt 
*'—  00 

=  (1/2t)  f  £,  (r,«)  Hy  (r,«) 

=  [/§  %/(8ir)l  J  *  I  F(«)  I  ^  (— )"  I  /(7.«)  I  ^  (14) 

'-00  cz 

where 
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|/(7.«)|2  =  (cz/u>)2  [U\  (7.«)  +  Ul  (7,M)]. 


(15a) 


which  is  suitable  for  p/a  =  s  >  \  and 

|/(7.«)|^  =  [1  +  Vl  (s,u)  +  V]  (5,«) 

-  2Vq  (s,u)  cos  {(7  +  u^/y)/2]  -  2  KiCj.k)  sin  {(7  +  M^/7)/2n,  (15b) 

which  is  suited  for  5  <  1. 

Except  for  5  =  0  and  j  =  1,  none  of  the  above  integrals,  (13)  and  (14),  can  be  evaluated  in 
closed  forms.  However,  for  small  values  of  s,  |  Iiy,u)  |  ^  given  by  (15b)  can  be  used  in  a  series  form 
in  powers  of  s.  On  the  otherhand,  for  large  values  of  s  Eq.  (ISa)  can  be  used  to  expressing  it  as  a 
series  in  powers  of  (1/s).  Similar  statement  holds  for  the  expression  of  the  electric  field  (13).  Let  us 
now  express  I(y,u)  for  various  values  of  s,  namely,  5  =  0,1,  ^  «  1  and  s  »1.  These  will  be 
found  to  be  very  useful  in  evaluating  £i(r,r)  and  P(r).  For  s  =  0,  it  is  easier  to  compute  I{y,u) 
from  (6b),  although  (7)  and  (10)  may  also  be  used.  Thus  for  f  =0  we  have 


2  1  I  •  2  1 

/(7.«)U -0  =  2(c2/«)  sm  exp  , 


|/(7.«)1,  =0  1^  =  4(cz/«)^  sin^  , 

4CZ 

V.  V 

For  s  «  1  using  (9),  (10),  and  (11)  one  finds 


u  —  ka^s/z. 

Note  that  for  very  small  s,  Jo(u)  —  —  0(s)  and  J2(u)  =  0(J^). 
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|/(7.«)l^  “  (z/k9  [l  +  4(«)  -  2  cos  (ka^/2z)JQ(u)  -  2j  sin  (ka^ /2z)  J i{u) 
+  ~  2Jo(u)J2(u)  +  2  cos  {ka} /2z)  J 2{u^  +  •• 


s  <  <  \. 


(17b) 


The  expressions  (17a)  and  (17b)  reduce  to  (16a)  and  (16b)  respectively  for  s  =  0.  For  ^  =  1 
either  (7)  or  (10)  may  be  used  and  the  results  are 


/(7.«),  -  1  =  7(7.7)  =  Hzc/ui)  exp  (17/2)  -  yo(7)]  /2.  (18a) 

and 

1 7(7.7)  I  ^  =  (cz  [1  +  4(7)  -  2  /o(7)  cos  i\/A.  (18b) 

For  very  large  values  of  s,  the  representations  (7)  and  (15a)  are  appropriate.  Thus  we  have  for 

j  >  >  1 


7(7,«)  “  (zc/«)  exp  (17/2)  [J\{u)/s  -  U2{u)/s^  +  0{\/s^)],  (19a) 

and 

|7(7,«)|^  =  {zc/w}^  j]{u)/s^,  s  »  1.  (19b) 

n.  BEHAVIOR  OF  FIELD  AND  ENERGY  DENSITY  IN  THE  FRESNEL  REGION  DUE  TO 
EXCITATION  BY  SINUSOIDAL  PULSES 

In  this  study  we  shall  consider  only  two  sinusoidal  current  pulses,  namely  single-cycle  sincdo' 
and  cosciiot.  where  ti>o  is  the  angular  carrier  frequency.  The  respective  spectra  of  these  pulses  are 

F,iu))  =  -  coo  (1  -  exp  (iti>T))/(w^  -  Wo).  (20a) 

for  the  sine  pulse  and 

Fc(u)  =  im(1  “  e*P  -  wo).  (20b) 

for  the  cosine  pulse,  where  T  »  Ix/uq.  As  the  frequency  u  increases  indefinitely,  one  finds 
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(21a) 


|>,(«)|  =0(l/a)2) 

and 


|/V(«)|  =  0(l/«).  (21b) 

A  A 

For  single-cycle  pulses  sintoo/  and  cosojo^,  F(w)  in  (13)  and  (14)  are  replaced  by  F,(w)  and 

A 

Fc(o)),  respectively. 

Since  we  wish  to  compare  the  energy  densities  of  the  radiated  fields  associated  with  these 
single-cycle  pulses  at  a  given  point  in  the  Fresnel  region  with  that  of  a  CW  sinusoidal  signal  (sinuof 
or  costdo/),  it  is  necessary  to  normalize  the  energy  content  of  these  signals.  Therefore,  we  shall 
assume  that  the  energy  of  each  of  the  single-cycle  sinusoidal  pulses  is  the  same  as  the  energy  per 
cycle  of  the  corresponding  CW  signal.  This  assumption  then  leads  to  the  following  expression  for  the 
square  of  the  magnimde  of  the  spectrum  of  a  CW  sinusoidal  (sina>or  or  cosojqO  signal. 

lFcw(t>t)  1  ^  =  (T^/(2tfo))  [3(«  -  a>o)  +  6  («  +  «o)3.  (22) 

where  S  is  the  Dirac’s  delta  function.  Replacing  |  F(w)  |  ^  in  (14)  by  |  Fcw(u)  | ,  the  energy  density  of 

the  sinusoidal  CW  excitation  can  be  expressed  as 


Pcw0  =  Po  («o/(2cz))2  1  Iiy,u)  U  =  ^  I  ^  (23) 

where  Pq  =  (/oVo/^Xt/uq)  is  a  constant,  which  has  a  dimension  of  Watts-sec  per  meter.  It  may  be 

noted  that  |  I(y,u)  |  ^  is  an  even  function  of  a. 


For  s  =  0,  using  (16b)  in  (23),  the  energy  density  along  the  axis  of  the  circular  disk  due  to  a 
CW  sinusoidal  excitation  is  given  by  the  well-known  expression  [3,5]. 


r  ^ 


Pcw(r)  1,-0=  Pcw(z)  =  Po  sin^ 


T 

8a 


where  a  -  z  /(4a^/\))  is  a  normalized  distance  along 
Next  from  (17b)  and  (23)  we  have  for  s  «  1 


(24) 

the  axis  of  the  circular  disk  and  Xq  —  2rc  /u>o. 
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Pcw(^  =  (/*o/4)  [1  +  4 


T5 

2a 


-  27(1 


TS 

2a 


cos 


T 

4a 


-  2sJx 


vs 

4a 


sin 


V 

4a 


+  s^ 


vs 

2a 


-27n 


vs 

2a 


vs 

2a 


+  27, 


vs 

2a 


cos 


r 

V 

4a 


+  -I. 


which  reduces  to  Eq.  (24)  when  5=0. 


For  5  =  1,  one  obtains  from  (18b)  and  (23) 


Pcw(r)  =  (7’o/16)[l  +  4 


2a 


-  2Jr 


r 

V 


2a 


cos 


r 

V 


2a 


Finally,  for  5  »  1  me  have  from  (19b)  and  (23) 


Pcw(P)  =  (Pq/^)  a 


vs 

2a 


/s^. 


(25) 


(26) 


(27) 


Before  we  proceed  to  compute  the  energy  density  integral  (14)  for  the  sinusoidal  pulses  let  us  now 
calculate  the  time  dependent  electric  field,  for  both  the  pulses  at  various  values  of  s.  For  this  purpose 
the  expression  (13)  will  be  evaluated  using  (7),  (10),  (16a)  and  (18a).  The  integration  technique 
needed  for  the  evaluation  of  (13)  is  discussed  in  the  Appendix  B.  Here  we  shall  simply  present  the 
pertinent  results.  On  the  axis  (5  =  0),  the  fields  Exs{z,t)  and  Exc{z,t),  excited  by  the  single-cycle 
sine  and  cosine  pulses,  respectively  are 


£«(z.O  =  ^0  [-  sin  Wot*  {C/(t*)  -  U(t*  -  T))J 


2  2  2 


Exc(z,t)  =  Eq  COSCJo^*  -  Uit*  -  7)) 


cos«o  (t*  -  ^)  mt*  -  -  V(t*  -  -  7)1 


(28a) 


(28b) 


where  Eq  =  /q  ijo/2  is  a  constant,  which  has  a  dimension  of  electric  field,  t*  =  t  -  z/c  and  U{f)  is 
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the  unit  step  function.  Note  that  (28b)  can  be  obtained  from  (28a)  by  differentiating  E„{z,t)  with 
respect  to  uot*.  This  should  be  expected,  since  E„  and  Ej,c  are  excited  by  single-cycle  sinuot  and 
costdQt,  respectively.  The  first  term  of  each  of  the  expressions  (28a)  and  (28b)  is  radiated  directly 
from  the  center  of  the  disk,  whereas  the  second  term  of  each  expression  is  radiated  from  all  points  on 
the  circumference  of  the  disk.  Since  the  point  of  observation  in  this  case  is  on  the  axis,  from  which 
any  point  on  the  circumference  is  equidistant,  the  radiated  field  from  each  circumferenical  point  is 
equal  and  all  of  them  add  nicely  to  a  single  term  (the  second  term  of  (28a)  and  (28b)).  However, 
when  the  observation  point  is  off-axis  (i.e.  s  is  finite  and  nom^ro),  its  distance  from  different  circum¬ 
ferential  points  is  different  and,  therefore,  an  infinite  number  of  terms  are  needed,  in  general,  to 
represent  the  total  off-axis  field  or  energy  density.  For  s  <  1  the  electric  fields  have  the  following 
representations. 


Exs(r,t)  =  -  £o  sm<oo/*  ((/(/*)  -  U(,t*  -  7)1 


+  £otE(-l)'"5^  /2m 

m  «0 


r  s 

rs 

2a 

j 


»(1  +  s^) , 

sin  «o  {^* - - -) 

4auo 


-  E  (-D^^^^'/im+l 

m  =0 


r 

IL 

2a 


,  ^  ir(l  -I-  „ 

cos  «o  It* - ;; - 11 

4a 


.  [[/{,*  _  5.0  +^)_  I  _  j;  j,#  _  _  7^]^ 

4acdo  4aa)o 


where  we  have  used  a^cjo/(cz)  =  ▼/(Za)  and 


(29a) 


Exc  (7,0-  -  Eq  cos  UQt*  I £/(/*)  -  U(t*  -  7)1 


+  Eq 


m  »0 


/2m 


COSOJQ  |t*  - 


y(l  +  5^)  I 
4awo 


m  ~0 


sin  cjQ  It* 


ZlLjLiil) 

4acao 


17 


.  £/{,*  _  _  t/{r*  -  -  71  .  (29b) 

4acao  4a(t)o 

In  this  case  tte  first  terms  of  (29a)  and  (29b)  also  correspond  to  the  radiation  directly  from  the  center 
of  the  disk.  The  rest  of  the  infinite  series  of  each  expression  represents  the  contributions  from  every 
point  on  the  circumference  of  the  disk.  For  s  =  0  each  of  the  infinite  series  reduces  to  a  single  term 
which  is  equal  to  the  respective  second  term  of  (28a)  and  (28b).  Here  again  (29b)  can  be  obtained  by 
differentiating  (29a)  with  respect  to  oiQt*.  For  5  =  1,  using  some  well  known  relations  involving 
Bessel  functions,  the  expression  (29a)  and  (29b)  can  be  reduced  to  (30a)  and  (30b)  respectively. 
Alternatively,  using  (18a)  into  (13)  together  with  (20a)  and  (20b)  the  relations  (30a)  and  (30b)  can  be 
obtained.  Thus  for  s  =  1  we  have 


EjcsiP  =  a,z,t)  =  (Eo/2)  [-  sin  <oot*{f/(r*)  -  U{t*  -  D) 


+  70  -  T^)  -  — )  -  -  —  -  71)  ,  (30a) 

2a  2a(i)o  ^*^0 


E,Ap  =  o,z,t)  =  (Eo/2)  [-  coswot*  \U(t*)  -  U(t*  -  71) 


+  7o  ^  cos  wo(t*  -  -:^)  U(t*  -  -^)  -  U(t*  -  -T)\  .  (30b) 

[^2aJ  2a<<7o  awo  “‘*’0  J 

Noting  that  t* - —  =  t  -  (z  +  2a^/z)/c  =  r  -  {(2a)^  +  for  2a /z  <  <  1,  one  may 

a«o 

attempt  to  point  out  the  origins  of  each  terms  of  (30a)  and  (30b)  in  the  following  manner.  Consider  a 
point  P  on  the  cylindrical  surface  of  radius  a,  the  base  of  this  cylindrical  surface  being  the  circular 
disk  radiator.  Let  ^  be  the  intersection  of  a  line  through  P  parallel  to  the  z-axis  and  the  circumfer¬ 
ence  of  the  circular  disk.  Then  PQ  is  the  z-coordinate  of  P.  Let  the  point  /?  be  on  the  circumference 
diametrically  opposite  to  Q,  i.e.  QR  =  2a.  Then  PR  =  ((2a )^  +  z^)'^^  =  z  +  2a^/z,  for 


2a  «  z.  The  first  term  of  (30a)  and  (30b)  shows  the  portion  of  the  field  at  P  originating  from  Q. 
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The  second  term  is  contributed  by  rays  from  all  other  circumferential  points,  the  point  R  having  the 
largest  distance  from  P. 

The  electric  fields  for  points  s  >  I  can  also  be  expressed  in  the  form  of  an  infinite  series  simi¬ 
lar  to  (29a)  and  (29b)  using  the  relations  (7)  in  (13).  However,  we  simply  present  the  dominant  term 
for  only  2/a  »  s  »  1. 

r  5  r  2 

r-  r  ,  F-  ^  ir(l  +  s^) ,  t(1  -I-  sY  , 

Exsir,t)  =  -  EQ(l/s)Ji  —  cos  Wo  (/* - ^ - )  U{t* - - - j 

2a  4awo  4awo 

V-  -/  L 

-  U{t*  -  -  71  (31a) 

4awo 

and 

r  1  2  r  2 

^  rs  ir(l  -t-  s^)  ^  t(1  -t-  s)  , 

Exc(r,t)  *  Eo(l/s)Ji  —  smwo  |t* - ;; - j  - r - ) 

2o  4awo  4a«o 

^  J  L 

-  Ult*  -  -  71  .  (31b) 

4a«o 

Let  us  now  evaluate  the  expression  for  the  energy  density  P{r)  given  by  (14)  for  those  special 
cases  used  in  electric  field  computations.  Before  proceeding  farther  let  us  reexpress  (14)  in  a  dimen¬ 
sionless  form  which  is  more  suitable  for  numerical  computation.  Then  for  both  the  sinusoidal  pulses 
we  have 

P,h  =  (Po/(2c^)l  5„'  (32» 

and 

PAr)  -  I<>o/(2a)'l  I I '  ‘‘i-  (32b) 

where 
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(33) 


H^,^,a,s)  =  r  -^0 


r 

2a 

V.  J 

exp 

4a 

c  j 

f</r- 


For  the  special  cases  where  closed  from  solutions  are  possible,  we  shall  use  (14)  together  with  (20a) 
and  (20b),  and  the  results  presented  in  Appendix  B  will  be  employed  without  further  explanations. 
Thus  along  the  axis  of  the  disk  where  s  =  0,  we  have  with  the  aid  of  (16b)  the  expressions  for  the 
respective  energy  densities. 


Ps(z)  =  Po  -  a) 


(1  “  (1  ~  -S— )  cos 

SOf 


T 

4a 


sin 


r  ^ 

T 

4a 


2t 


-)t/(a  -  1/8) 


and 


Pciz)  =  Pq  [f/(l/8  -  a) 


+  “  (1  “  IT-)  cos 

8a 


r 

■K 

4a 

V..  J 


sm 


4a 

-^^)C/(a  -  1/8) 


(34a) 


(34b) 


It  follows  from  (34a)  and  (34b)  that  for  large  values  of  a,  for  which  ir/4a  <  <  1,  P^{z)  and  P^(z) 
behave  respectively  as 


f’,(z)~  Pq  ir^/(32a^). 


(35a) 


and 


Pr(z)-  ^0/(4a).  (35b) 

Relations  (35a)  and  (35b)  show  that  in  the  far  end  of  the  Fresnel  region  the  energy  along  the 
axis  due  to  a  single-cycle  sine  pulse  decreases  like  1/z^,  whereas  the  single-cycle  cosine  pulse  causes 
the  energy  to  decay  as  \/z,  which  is  much  slower  than  that  for  the  sine  pulse.  The  corresponding 
CW  energy  density  (see  (24))  behaves  like  1/z^,  which  follows  from 
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Pcwiz)  ~  Po  ir2/(64a2).  (35c) 

Although  both  P,  and  Pew  dn^  like  1/z^  in  the  far  end  of  the  Fresnel  region,  the  energy  density 

associated  with  the  single-cycle  sine  pulse  is  double  that  of  CW.  This  shows  that  both  the  single¬ 
cycle  sine  and  cosine  pulse  can  really  enhance  the  Fresnel  region  when  conqnred  with  the  CW  signal. 
Noting  the  high  frequency  behaviors  of  the  spectra  of  the  sine  and  cosine  pulses  (see  (21a)  and  (21b)), 
it  can  be  stated  that  the  higher  the  high  frequency  content  of  a  pulse,  the  larger  is  the  extension  of  the 
Fresnel  region.  Although  a  single-cycle  cosine  pulse  is  unrealizable  (similar  to  an  ideal  rectangular 
pulse),  the  result  illustrates  the  fact  that  the  shorter  the  rise  time  (which  corresponds  to  amount  of 
high  frequency  content)  of  a  pulse,  the  greater  is  die  distance  along  which  the  energy  density  remains 
highly  collimated.  This  observation  suggests  that  if  some  realizable  pulses  can  be  created  having  the 
frequency  spectra,  the  behavior  of  which  lies  between  l/o)^  and  l/w,  for  very  large  u,  such  pulses 
are  potential  candidates  for  the  enhancement  of  the  Fresnel  region. 


Next  we  present  the  off-axis  energy  density  in  the  vicinity  of  the  axis  (i.e.  for  s  «  1),  using 
(17b)  into  (14),  in  the  following  manner  [see  the  (terivations  of  the  integrals  (2)-(4)  in  Appendix  B] 

16 

Ps(P,z)  =Po  L  fs.i,  s  <  1/2  (36a) 

i-l 

16 

P^(p,2)=  Pe  Z  5  <  1/2  (36b) 

i-l 

where 


r 

Is,  1  *  -(1/4)  n(2/T)sin 

T 

4a 

V-  j 

—  (l/2a)  cos  j 

IT 

4a 

L  J 

Uo 

ITS 

2a 

+  (s  /a)  sin 


]U(a), 


r 

r 

It.  2  *  (J/4)n(4/»)  cos 

y 

4a 

^  J 

+  (l/2a)  sin 

JL 

4a 

xs 

2a 

^  J 

(37a) 


(37b) 
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r 

4,o=(jV2)y?  t/(a-5/2). 

w  ^ 


^  ^2  ^  Vict-sH), 


U^/8)n-(6/T)sm  ^  +2(2--^)  cos  ^ 


+  (j/a)y,  ^  sin  ]C/(a  -  1/8  -  s/4), 

2a  4a 


(f.  13  =  (^^/8)  [|(6/t)  sin 


(37m) 


-(j/of)^!  sin  ]i/(a), 

2«r  4a 


/M4»(l/8)n2/T)sin  +2(2-;^)  cos  —  l/o  2^r 


+  (s/a)y,  ^  sin  ^  ]J/(l/8  -  J/4-2), 

2a  4a 


15  =  /8)ll  -  (4/t)  cos  —  +  2(2  - 


1  .  .  IT  ,  ,  « 


+  (j/a)yo  ]£/(l/8--j  -  a), 

2a  4a  4 


16  =  -  (^V8)  [((6/ir)  sin  —  +  2(2  - 


1  .  »  1  I 

— )cos  —  1^2  T" 
4a  4a  2a 


(5/a)y,  ^  sin  ^  1  I/(l/8  -  5/4  -  a). 

2a  4a 


4.1  =  (l/4)l{(2/ir)sin  +(l/2a)cos  1 -/o 
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-  (j/o)/) 


}U{a) 


+  (s  /a)  Jq 


+  is /a)  J\ 


—  is  /2a)  J\ 


—  is  /2a)  J 1 


+  is /2a)  Jq 


ITS 

2a 


sm 


r 

T 

4a 

j 


Ic.2  =  -  is/4)[il/2a)Ji 


xs 

2a 


sin 


r 

T 

4a 


(38b) 


xs 

2a 

J 


cos 


T 

4a 


]!/(«). 


Ic.3  =  -  is^/4)  [(-  (2/t)  sin 


r  "N 

X 

4a 


+  (l/2a)  cos 


r 

X 

4a 


xs 

2a 


(38c) 


xs 

2a 


sm 


X 

4a 


]f/(o). 


4.4  =  (l/2)£/(a), 


4.5  =  i\/2)Jl 


xs 

2a 

J 


Via), 


(38d) 

(38e) 


4.6  =  (1/4)  I((l/ir)  sin 


r 

ir 

4a 


—  2(2  —  l/4a)  cos 


r  ^ 

X 

4a 


Uo 


xs 


(380 


xs 

2a 


sm 


4a 

V.  ^ 


]Uia  -J  -  s/4). 


4,7  =  -  (1/4)  [{!/»)  sin 


r 

n 

4a 

V.  ^ 


+  (2  +  l/4a)  cos 


r  'N 

T 

4a 

Vw  ^ 


Uo 


r  ^ 


2a 


(38g) 


xs_ 

2a 

J 


sm 


4a 

J 


1  Via), 


4.8  =  -  is/4)li2  -  l/4a)/, 


xs 

2a 

V.  y 


sm 


r 

X 

4a 


(38h) 


]  Via  -  1/8  -  5/4), 
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-  (s  /2a)  Jo 


—  (s  /2a)  J\ 


+  {s  /2a)  J  \ 


+  (i  /2a)  J 1 


+  (s/2a)  Jq 


4.9  =  -  (J/4)  [(2  +  l/4a)  J, 


TS 

2a 


sin 


jr_ 

4a 


(38i) 


■KS 

2a 

V.  ^ 


cos 


4a 


]!/(«), 


4, 10  =  (s^n)  J\ 


rs 

2a 

V. 


Via  -  s/2). 


(38j) 


4,11  —  s^ 


r 

vs 

2a 


Jl 


r 

vs 

2a 


U  {a  —  s  /2), 


(38k) 


4,12  =  (^^/4)  [{(l/x)  sin 


+  (2 - ^)cos 

^  -N 

TT 

)-/2 

^  "N 

w 

4a 

v-  y 

4a 

r 

4>> 

P 

2a 

(381) 


r 

vs 


2a 

J 


sin 


r  ^ 

V 

4a 


]  Via  -  1/8  -  ^). 

4 


4,13  =  (5^/4)  [(-  (l/x)  sin 


X 

+  (2  +  -^)  cos 

X 

)-/2 

JEl 

4a 

4a 

4a 

L  J 

2a 

L  J 

(38in) 


r 

XJ 

2a 

V.  y 


Sin 


r 

V 

4a 

y 


]U{a), 


4.14  =  (1/4)  [|-  (l/x)  sin 


r 

X 

4a 


+  (2  -  — )  cos 
4a 


4a 


\Jo 


vs 

2a 

V  y 


(38n) 


r 

vs 

2a 

y 


Sin 


r 

V 

4a 

V.  y 


'c.l5  -  -  (»/4)(C  -  l/4<»)7| 


r 

VS 

2a 

y 


Sin 


r 

X 

4a 

v  y 


(38o) 


vs 

2a 

V.  ^ 


cos 


4a 

V. 


1  17(1/8  -  y  -  a), 
4 


4. 16  =  -  (•*^/4)  [{(l/x)  sin 


X 

4a 


r  ^  r  ^ 


+  (2  —  l/4a)  cos 


X 

4a 

v.  y 


\h 


VS 

2a 

V.  y 


(38p) 
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-  (5 /2a)  J\ 


r  ~\ 

TS 

4a 

V  y 


sm 


4a 

v.  y 


]  U(l/&  -  S/4  -  a). 


Note  that  for  j  =  0,  the  relations  (36a)  and  (36b)  reduce  to  (34a)  and  (34b)  respectively.  Then 
for  s  =  1,  using  (18b)  in  (14)  we  have 


P,(p  =  a,z)  =  (Po/8)  tl  +  4 


2a 


-  2Jo 


r 

r 

2a 


cos 


r  ~\ 

T 

2a 

y 


+  (l/2a)  (Jo 


r 

K 


2a 

V  y 


cos 


''i'' 

2a 

V.  y 


-Jx 


r 

V 

2a 

V-  y 


sm 


r 

T 


2a 


-  ir 


r  ~\ 

ir 

2a 


sui 


r  ^ 
jr_ 

2a 

V..  y 


/tJ,  a  S  1/2. 


and 


(39a) 


P,(p  =  a,z)  =  (Po/8)  [1  +  4 


2a 

v  y 


-  2Jo 


2a 

y 


cos 


2a 

v.  y 


(l/2a)  Uo 


r 

JL 

2a 

v.  y 


COS 


r 

y 

2a 

V  y 


-h 


SUI 


2a 

V.  y 


+  Jq 


r  ~\ 

T 

2a 


Sin 


r 

y 

2a 

V.  > 


/y],  a  2  1/2. 


(39b) 


Finally  for  s  » 1,  use  of  (19b)  in  (14)  yields  the  following  expressions  for  the  energy  densities. 


P,  (p,z)  =  Po  4 


KS 

2a 


/(2s2)  =  P,  (p,z),  \«5  <  2a. 


(40) 


The  conditions,  such  as  1  «  5  <  2a,  f/(a  -  1/8  -  i/4)  etc.,  which  appear  with  various 
expressions  indicate  the  region  of  validity  of  that  expression  in  question.  These  conditions  are  noth¬ 
ing  but  the  requirements  of  the  convergence  of  certain  integrals,  which  had  to  be  evaluated  in  obtain¬ 
ing  that  expression.  In  addition,  it  is  also  necessary  to  remember  that  the  validity  of  the  approxima¬ 
tion  in  (2a)  used  for  Fresnel  region  requires  a^/(2z)  «  z,  ap/z  «  z  and  p^/(2z)  «  z,  which 
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are  equivalent  to  4a(a  /Xo)  »  ^  and  4a(q  /Xo)  >  >  s  /^,  for  all  values  of  s  and  a.  These  later 
conditions,  must  therefore,  also  be  imposed  simultaneously  with  those  mentioned  above. 
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APPENDIX  B 


EVALUATION  OF  CERTAIN  INTEGRALS 

Special  cases  of  the  following  integrals  involving  Bessel,  exponential  and  trigonometric  functions 
appear  in  the  Appendix  A. 


II 

F/(u)  J„(/3<j)  e  ~  du. 

(1) 

OD 

=  J.  . 

1  Fi  (w)  1  ^  Jn  03<o)  Jm  (/3«)  dw. 

(2) 

>.2,/  =  J  ' 

OD 

1  Fi  (u)  1  ^  J„  (|3u)  cos  Oo  w)  du. 

(3) 

00 

-<3,,  =  5. . 

1  F/(w)  1  ^  J„  (fiu)  sin  (00  w)  dw, 

(4) 

where  and  o  are  positive;  and  I  =  s  for  the  single-cycle  sine  pulse  and  /  =  c  for  the  single¬ 
cycle  cosine  pulse.  The  frequency  spectra  of  the  single-cycle  sinusoids  are  given  by 

F,(o))  =  2/  «o  sin  (uT/2)/(<a^  -  w§),  (5) 

Ff(«)  =  — 2«  sin  (cjT /2)/(«^  —  a>o),  (6) 

A  0S 

where  uqT  =  2ir.  Since  the  hinctions  F,(w)  and  F^fw)  are  bounded  (or  finite)  at  a>  =  ±  these 
points  along  the  real  axis  of  w  are  not  the  singular  points  (or  poles)  of  any  of  the  above  integrals. 
However,  the  factor  sin(ur/2)  has  different  exponential  behavior  in  the  upper  and  lower 

halves  of  the  complex  (i>-plane.  Therefore,  it  is  necessary  to  split  this  factor  into  exponential  form, 
(e'"^  -  l)/(2i),  and  then  consider  each  term  of  these  factors  separately  for  evaluating  the  resulting 
integrals.  However,  this  procedure  introduces  pole  singularities  at  u  =  ±  coq  on  the  real  (>)-axis  in 
each  of  the  separate  integrals.  In  order  to  avoid  this  difficulty,  we  replace  the  integration  path  from 
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—  00  to  +00  in  the  original  integrals  (1)  to  (4)  by  a  path  (or  contour)  C  (Fig.  2),  which  also  runs 
along  the  real  axis  except  at  w  =  4:cjo,  where  upward  (i.e.  above  the  real  axis)  indentations  are 
made,  thus  avoiding  the  path  running  through  <>>  =  +  uq.  This  is  permissible,  since  the  original 
integrals  before  spliting  the  functions  F/(<i>)  or  |  Fiiu)  |  ^  into  exponential  forms,  are  analytic  in  the 
neighborhood  of  w  =  +  coq,  an  appropriate  deformation  of  the  p)ath  of  integration  around  these  points 
will  not  change  the  value  of  the  integrals.  Therefore,  the  original  path  along  the  real  axis  for  all  the 
integrals  (1)  to  (4)  is  replaced  by  the  contour  C,  described  above. 

Before  evaluating  (1)  to  (4)  it  will  be  very  helpful  to  consider  the  following  auxiliary  integrals 
involving  Bessel  and  Hankel  functions  of  integer  orders. 


-  1  — 5 - ff  s  jS  >  0, 

••c  ((0^  -  4jg) 


(7) 


^  =  f  «  — ^ - r-  exp  da,  a  ^  0  >  0, 

c  {uT  -  Wo) 


(8) 


L  exp  i-iia)  d«,  8  >  /3  >  0, 


(9) 


J„0O}) 


=  \  a^  — - r-r-  exp  (-i8w)  da,  5  i  jS  >  0, 

(u,2  -  u§)2 


(10) 


^  t  03a>)  da,  0  >  0,  0  <  m  £  n, 

{up-  -  air 


(11) 


=  f  0>O,  Osmsn+2, 

C  (w  —  Uq) 


(12) 


J„(0a) 


=  1  — 5 - FT  Jm0<^)  e*P  i-ioa),  a  2:  2/3  >  0, 

Jc  {u?  -  air 


(13) 


Jg  =  [  a^  .  — 2  2  Jm  03«)  exp  i-iaa),  a  >  2  0  >  Q. 

C  (or  —  uq) 


(14) 
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The  contour  C  of  the  above  integrals  (7)-(14)  can  be  joined  by  an  infinite  semi-circle  in  the 
lower  half  of  the  complex  deplane.  The  integrands  on  the  infinite  semi-circle  vanish  either  exponen¬ 
tially  or  at  least  as  The  branch  cut  of  (0u)  in  (11)  and  (12)  is  chosen  along  a  ray  (start¬ 

ing  at  the  origin  u  =  0)  in  the  2nd  guadrant  of  the  to-plane,  making  a  very  small  angle  with  the  nega¬ 
tive  real  axis,  so  that  the  branch  cut  does  not  touch  the  indented  contour  C.  In  order  to  avoid  the 
branch  cut,  the  contour  C  in  (11)  and  (12)  is  indented  also  below  the  origin,  although  they  are  not 
shown  in  Fig.  2.  It  may  also  be  noted  that  imposition  of  various  convergence  conditions,  such  as 
osm<n,  2,  a>/3>0,  5&j8>0  and  a  s  2  j8  >  0,  the  behaviors  of  the 

Bessel  and  Hankel  functions  in  the  complex  oj-plane  played  a  role.  With  these  preliminary  observa¬ 
tions,  the  preceeding  integrals  can  be  evaluated  easily  by  the  application  of  the  calculus  of  residue.  In 
particular,  if  N((i})  represents  any  of  the  numerators  of  the  integrands  of  to  then  the  value  of 
those  integrals  can  be  expressed  as 


(-2tj) 


du 


Niw)/(<a  —  u>or 


^  +  —  ■|N(w)/(a)  -f-  (Oof  i- 

<tf=—  Uq  L  J  U  =  U[) 


(15) 


Following  the  procedure  outlined  above,  the  integrals  (7)  to  (14)  can  be  evaluated  as  presented 
below. 


r(-2»/«o)  Jn  (/3«o)  sin  (wq  o)  U  (a  -  0), 
I  for  n  =  0,  or  even  integer. 


(16a) 


[-(2x1  /«o)  J„  (pcjo)  cos  («o  o)  U(a 
for  n  =  odd  integer. 


(16b) 


where  U{x)  — 


1,  forx  s  0 
0,  forx  <  O’ 


(17) 


r(-2xi)  (/S  «o)  cos  (wo  <r)  U{a  -  /3) 

I  foni  =  0,  or  even  integer. 


(18a) 
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'(-2t)  J„  (fiwd)  sin  (uq  a)  U(o  -  0) 
for  n  =  odd  integer. 


(18b) 


(»/«o)  [(«  +  1)  03  Mo)  sin  (6  ti)o)  -  /3  «o  •^n-i  05«o)  sin  (&i)o) 

^  =  «  -5mo  y„  (i8mo)  cos  (5mo)1  U  (8  -  0),  (19a) 

for  n  =  0,  or  even  integer, 

r  . 

(it/mo)  K"  +  1)  Jn  (0o>q)  cos  (S  Mo)  -  0Uo  Jn-l  (/5mo)  cos  (6mo) 

=  -  +  «  MO  03mo)  sin  (8  mo)]  U(8  -  0)  (19b) 

for  R  =  odd  integer. 

(t/mq)  [(n  -  1)  Jn  03mo)  sin  (Smq)  -  0u>o  J„-i  {0u>o)  sin  (Smq) 

^  =  -  -  8o}QJn  03mo)  cos  (5mo)1  U  (8  -  0)  (20a) 

for  R  =  0,  or  even  integer, 

(iT/to^)  [(r  -  1)  J„  (0o3o)  cos  (5mo)  -  0UO  J„.i(0u)o)  cos  (6  mo) 

=  •  +  5mo  y„  03mo)  sin  (Smq)]  U  (8  -  0),  (20b) 

for  R  =  odd  integer. 

(iir/ul)  ({(n  +  m  +  l)y„  (/Smo)  -  0u)o  Jn-\  (/3«o))  Jm  (/S«o) 

^  —  0U0  J„  (0uq)  J„~i  (jSmo)]  (21a) 

for  R  +  Ri  =  0,  or  even  integer, 

(r/wo)  U(n  +  m  +  1)  y„  (0wo)  -  0uo  •/n-i(/3“o))  OSmq) 

=  -  —  0(t)QJ„  (/3mo)A/^_i  (^mq)]  (2lb) 

for  R  +  RI  =  odd  integer, 

where  Nn(0UQ)  and  N„-i{0(>)o)  are  Neumann’s  functions. 

Taking  now  the  real  and  imaginary  parts  of  (11)  and  (21)  we  have 


J„(0u)  J„  (0w)dw 

-  Ic  ^y‘ 

{b,for  R  +  RI  =  0,  or  even  integer, 
0  £  RI  s  R, 


(22a) 
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(t/«o)  [|(«  +  ffl  +  1)  y„  (/3<Jo)  -  /3<i)o  Jn  --t  (/3f>o)Wm  03«o) 
-  iScoo  (l8«o)  ^M-1  03«o)l 
for  «  +  m  =  odd  integer, 

0  s  m  :s  n. 


^5b  =  [ 


,  J„  (0ui)  N„  (/3<j)  du) 

c  (td^  _  4j2)2 


(t/(ijo)[(“  (n  +  m  +  1)  J„  (/Stoo)  +  |8ci>o^n-i  (/3<«>o))^m  (/3uo) 

+  0UqJ„  0<do)/m_l  (j3ci>o)l 
for  n  +  m  =  0,  or  even  integer, 

0  £  m  <  R, 


=  0,for  R  +  Ri  =  odd  integer.  0  <  ri  <  r. 


(-iT/cdo)l{-  (R  +  RI  -  1)  J„  (jStdo)  +  jSwo  Jn-l  i^Uo)]  Jm  (/3wo) 
+  0<i>o  J„  (Pwq)  Jm-\  0<«>o)l* 
for  R  +  RI  «  0,  or  even  integer 
0  S  RI  s  R  +  2 


(-T/uo)  [(-  (r  +  RI  -  l)y„  (jSuo)  +  /3wo  ^-1  (|3«o))  (/3u)o) 

=  -  +  /3<|J0  /n  (/3«o)  A'm-l  (/3uo)l 

for  R  +  RI  =  odd  integer,  0  S  ri  s  r  +  2. 

Equating  real  and  imaginary  parts  of  (12)  and  (24)  we  have 


w^y„  (M  /„Oc3)  ^  ^ 

■  *c  <<«.  0  s  m  s  n  +  2 


(fiT  -  «o) 


=  ■<  0,  for  R  +  RI  =  0,  or  even  integer 


(  — ir/<«)o)[(“  («  +  W*  ~  i)«^n  (^<^)  +  "^n-l  (/3wo))^m  (^<*10) 

=  '  +  /SciJo  (/3wo)  - 1  (^t^)] 

for  R  +  RI  =  odd  integer. 
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‘^  =  J 


.  <J?  J„  (0ui)  N„  (fiu)du) 


,  0  £  m  S  n  +  2 


(t/«o)[-  {n  +  m  -  l)J„  (/3«o)  +  0uioJn-i  03«o))  Jm  Owo) 
+  fiuQ  J„  03uo)  im-i  03«o)l 
for  «  +  m  =  0,  or  even  integer. 


=  -<  0,  for  n  +  #fi  *=  odd  integer 


(t/wo)  [(«+«  +  iy„  (/Scjo)  Jm  Owo)  sin  (ecoo) 

-  ^<1>0  {^n-l  (0<>>o)  Jm  (/3<Jo)  +  -Ai  (/3ci)o)  Jm-\  O^OJo))  Sin  (atiJo) 

-  ao)Q  J„  (/3«o)  •/«  (/3wo)  cos  (fftdo)l  U{a  -  2/S), 
for  n  +  m  =  0,  or  even  integer. 


(ir/wo)  [(«+»!  +  1)  y„  (jScJo)  •/m  (/3wo)  cos  (a  wq) 

-  l-fn-l  i$t^)  •fm(/3‘<>o)  +  "fw  03a>o)  ^m-l  (i8«o))  COS  (ff  (Jq) 
+  <r  coq  y„  03«o)  03«o)  sin  (a  o^)]  £/  (e  -  2/3) 

for  n  +  m  ==  odd  integer. 


(-  iT/«o)  n-  (n  +  m  -  •!)/„  O«o)  +  /3«o  O«o)) 

~  (^<*>o)  Jm-l  (/3<^)1  U  (a  —  2/3) 

for  n  +  m  =0,  or  even  integer. 


(-t/«o)  [|-  (/I  +  ffi  -  1)  (/3uo)  +  /3«o  4-1  (/8«o))  A^m  (/8a>o) 

=  '  +  /3<jo  ^  (/3«o)  A^m-i  (/3«o)]  £/  (<r  -  2/3) 
for  /I  +  m  =  odd  integer. 


In  addition,  the  following  relations  can  be  easily  established. 


.  ,  J„  (/3w)e“^d« 


Ci)  J„(0u)  e"^d(j) 
-  Wo) 


=  0,  <7  >  /3  >  0 
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(31) 


,  ,  y„OSa>) 

(.^-4)^  =0.i>3>0 


^2  =  1 


(jS^  J„(0u)  e'^da 
C  (a>2  - 


=  0,8  >  iS  >  0, 


(32) 


“  *c  — - 0,s>0,0sm=n 

.  V.:2,2 - =0.^>0,0^.s„+2. 


(33) 


(34) 


For  (33)  and  (34)  the  contour  C  is  also  indented  above  the  origin  and  the  branch  cut  of  is 

chosen  along  a  ray  from  u  =  0  to  —  oo  slighdy  below  C  in  the  3rd  quadrant. 


- o.<,>2«>o 


(35a) 


Ai  -  L 


7„(^«)  y„(Me'"“  du 


=  0,  ff  >  2|8  >  0. 


(35b) 


'C  («2  -  «g)2 

For  all  of  the  integrals  (29)  to  (35),  the  contour  C  can  be  closed  by  an  infinite  semi-circle  in  the 
upper  half  of  the  complex  t^-plane.  Since  all  of  the  integrands  vanish  on  this  infinite  circle  and  the 
integrands  are  analytic  inside  the  closed  contour,  the  integrals  vanish.  Furthermore,  we  have 


•^17 

•^18 

J^g 


*^0 


•^1 


r  _ ^ —  =  0 

(0,2  -  a,§)2 


{  — 
ir  /.  .2 


'c  -  a,§)2 


=  o,r  >  0 


i, 


,  -iwr 


c  (0,2  -  o,g)2 


=  -  vT/uio,T  >  0 


f  — i 

ir  t.  J- 


—iut 


r  '\ 


C  (o,^  -  0,6) 


ddl  =  — 


2ir 

«o 


sin  0,0  ^  '  t^(^) 


I, 


fa,  e 


-itH 


'c  (fa,2  _ 


=  -  2xi  cos  Wot  U  (t) 


(36) 


(37) 


(38) 


(39) 

(40) 


35 


In  view  of  (7)  and  (35),  we  have  the  following  relations. 


g  t  Jn  (iSt*))  JmWa)  COS  {au))d<a  .  ^ 

- ^57747 - J,n.,>2e>0. 


-  ,  J„  (0<a)  J„  (003)  sin  (a  «)  da 

- j,,m.o>2s>0. 

Similarly,  from  (9)  and  (31),  we  have 


.  .  J„(0O))  COS(5  03)d03  ^  ^  , 

■■  *c  -  4)2 - JS/2.  a  >  ^  >  0 

.  f  J„  (0a)  sin  (6  «)  d  « 

*^25  =  L - 71 - 2I -  =  -  ‘^3/(2i).  6  >  /3  >  0. 

C  (a-‘  -  Wo) 

From  (10  and  (32),  we  have 


=  f 


a^  J„(0a)  cos  (b  a)  d  a 


c  (w^  -  wo)^ 


=  Ji/2,  6  >  jS  >  0, 


-  f  0^  J„  (0a)  sin  (b  a)  d  a  , 

“  *c - ^4^455 - -  -  ■</ao.  a  >  «  >  0, 

Similarly,  from  (14)  and  (35),  the  following  relations  can  be  obtained. 


As  -  i 


a^  J„(0a)  J„(0a)  cos  (a  a)  d  a 


(a^  -  wo)^ 


=  Ji/l,  o  >  20  >  0, 


A9  -  L 


a^  J„  (0a)  Jm(0a)  sin  (a  aa)  d  a 
(a^  -  wo)^ 


=  —  ^f%/(2i),  o  >  2a  >  0. 


Using  now  the  preceeding  results  we  present  the  integrals  (1)  to  (4)  for  /  =  s  and  c  in  the  following 


manner. 


Jf,s  =  [“  WJn  (/3w)e-'"“dw 
^  -  00 


r2T  j„  (0  0 
~  i  for  n  =  0 


Wo)  sin  (wo  a)  [£/  (a  -  jS)  -  C/  (ff  -  |8  -  T)] 


or  even  integer. 
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Ivi  J„  08  Wo)  cos  (wo  a)  lU{a  -  jS)  -  t/  (a  -  /3  -  D] 
for  n  =  odd  integer. 


(49b) 


=  f  Fc  (w)y„  OSw)^-"^^/ w 

•  — OD 


f- 


T  y„  03  Wo)  cos  (<r  wo)  {Via  -  &)  -  Via  -  &  -  T)\ 
for  n  =  0,  or  even  integer. 


'■-2x1  J„  03  Wo)  sin  {a  wo)  {Via  -  /3)  -  £/  (o  -  /3  -  7)] 
for  n  =  odd  integer. 


(50a) 

(50b) 


.J  =  J  I  I  ^  Jn  03w)  J„  03w)  w 

r(2x^/wo)  Jn  (/3wo)  y„  03  Wo)  1/  (x/wo  -  /3) 
I  for  n  +  m  =  0,  or  even  integer. 


(51a) 


r 


-(2x/wo)  [{-(n  +  m  +  1)  y„  03  Wo)  +  /3wo  y„_i  03wo))A(„  03  wo) 

+  ^  wo  y„  03  wo)  )V„_j  (/3  wo)lf/03) 

+(ix/wo)  [-  in  +m  +1)  J„i0uo)  J„  03wo)  +  jSwo  y„_t  (|8wo)  O?  t^) 
+  /3  Wo  y„  (/3  Wo)  y„_i  03  wo)l  l/(x/wo  -  |8) 
for  n  +  m  =  odd  integer. 


(51b) 


tei,c  -  J  I  ^c(«)  I  ^  y^  03  w)  y„  (/3  w)  </  w 

■(2x^/wo)  y„  (j8  Wo)  y^  03  Wo)  c/  (x/wo  “  iS) 
for  fi  +  m  =  0,  for  even  integer 


(52a) 


2x/wo)  [((n  +  m  +  1)  y„  03  wo)  -  /3wo  J„-\  (jS  wq))  N„  03  wo) 

-  /3  Wo  y„03wo)  lV„_i  (0wo)]  t/  03) 

*  -  (ix/wo)  Il(n  +  m  -  1)  y„  03  Wo)  -  /S  Wo  y„  _1  (/3  wo))  y„  (/3  wo)  (52b) 

-  ^  Wo  y„  (/3  wo)  y,„_i  03  wo))J  V  (x/wo  -  /3) 
for  n  +  m  —  odd  integer. 


'Ki.s  =  J  I  ^5(“)  I  ^  y^  03  w)  cos  OSo  w)  w 

•  —  CD 

=s  (x/2wo)  ((n  +  1)  y„  (jSwo)  sin  03o  wo)  —  P  uq  J„^\  (/3wo)  sin  ifio  uq) 


+  (2t  -  ^  Wo)  y„  03wo)  cos  (/3o  wo)l  U  (2ir/wo  -  ~  0) 

+  (»/wo)  l(n  +  1)  /„  03  Wo)  sin  03o  <^)  -  jS  wq  _  i  (/3  wo)  sin  {0q  uq) 

~  00  <t)Q  Jh  (ff  «o)  cos  (^0  ‘*>o)]  Ui0Q  —  0) 

+  (t/2wo)  («  +  1)  03  Wo)  SfO  (/3o  Wo)  +  j3  Wo  y„  _i  08  «^)  sin  (0o  wo) 

+  00  4^“  2*-)  y„  08  «o)  cos  03o  t^o)!  U  (0q  -  2t/wo  -  0) 

+  (t/2wo)  [-  in  +  1)  J„  03wo)  sin  i0Q  wq)  +  jS  wq  y„_i  (0  uq)  sin  (0o  wq) 

+  (2t  +  i3owo)  J„  03  Wo)  cos  03o  <^o)]  U  (2t/wo  +  0o  -  0) 

when  rt  *  0,  or  even  integer.  (53a) 

-  (ix/2wo)  [-  (/I  +  1)  y„  03  Wo)  cos  03o  ci^)  +  /3  wq  y„  _  i  03  wq)  cos  (00  wq) 

+  {2t  -  jSowo)  J„  08  «o)  sin  (0o  wo)]t/  (2x/wo  -  0o  -  0) 

+  (Jx/wo)  t(/i  +  1)  03  Wo)  cos  03o  tJo)  “  /8  Wo  y„  -1  (0  wq)  cos  (/3o  wo) 

+  00  Wo  y„  (/3  Wo)  sin  (0o  wq)]  U  (0o  -  0) 

+  (jx/2wo)[-  (n  +  1)  J„  (/3wo)  cos  03o  wo)  +  j3  wq  y„  _  j  03  wq)  cos  (0o  wq) 

+  (/3o  Wo  -  2x)  J„  03  Wo)  sin  (0o  wq)]  U  (0o  -  2x/wo  -  0) 

+  (ix/2wo)  [-  (n  +  1)  y„  03  Wo)  cos  (0o  wq)  +  0o  uoJ„-t  (8wo)  cos  (0o  wq) 

-  (2x  +  00  Wo)  y„  (/3  Wo)  sin  03o  wo)I  U  (2x/wo  -  0o  -  0),  (53b) 

when  n  =  odd  integer 

‘^2.c  =  I  I  ^f(")  I  ^  03«)  cos  03o  wo)  y  w 

— 00 
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—  (t/2mo)  [(«  ~  1)  /n  03  «o)  sin  03o  <«)o)  -  /S«o  Jh-\  08uo)  sin  OSq  «o) 


+  (2t  -  /So  «o)  05  «o)  cos  OSo  too)]  (2T/too  -  0o  ~  0) 

+  (t/wo)  [(n  —  1)  J„  OStoo)  sin  /3o<oo  ~  /3  too  -^n-i  08  too)  sin  (/So  too) 

-  /3o  too  08  «o)  cos  03o  too)]  03o  -  |3) 

+  (T/2too)  [-  (n  -  1)J„  03  too)  sin  (^0  too)  +  S  too  •/»-!  08  too)  sin  (/So  too) 

+  (^0  t^  -  2t)  y„  (S  too)  cos  (So  too)]  U  (So  -  2T/too  -  S) 

+  (T/2too)  [-  (/i  -1)  J„  (Stoo)  sin  (So  too)  +  i8  «o  y„_i  (S  too)  sin  (So  too) 

+  (2t  +  So  too)  -/n  (/8o  too)  cos  (So  too)]  (2ir/too  +  So  -  /8) 
when  «  =  0  or  even  integer, 

*  (jt/2<oo)  ]-  (n  -1)  J„  (S  too)  cos  (So  too)  +  S  too  ^-i  (0  too)  cos  (So  too) 

+  (2t  -  So  too)  A  03  too)  sin  (So  too)]  U  (2x/<oo  -  /3o  -  S) 

+  (ix/too)  [(n  -l)/„  (S  too)  cos  (So  t^)  -  |8  Wo  (/3too)  cos  (So  too) 

+  00  “0  -^n  (^  too)  sin  (So  too)]  (^0  ~  0) 

+  (/x/2coo)  [-  (n-1)  J„  (Stoo)  cos  (So  too)  +  Stoo  yn-i  (S  too)  cos  (So  too) 

+  (^0  "0  -  2x)  y„  (Stoo)  sin  (So  too)]  (/3o  -  2x/too  -  0) 

+  (jx/2<oo)I-  (/I  -1)  y„  (Stoo)  cos  (So  too)  +  Stoo  y„-i  (Stoo)  cos  (So  too) 

-  (2x  +  So  too)  y„  (S  too)  sin  (So  too)]  U  (2x/<oo  +  So  -  /3), 

when  n  —  odd  integer. 


(54a) 


(54b) 
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•^3.i  =  J  I  ^j(")  I  ^  03  u)  sin  (/So  «)  rf  t.) 

=  (it/2<ijo)[-  (n  +  l)/„  03t*)o)  sin  (/3o  «o)  +  /5  -/n-i  (iSwo)  S“  O^o  ‘^) 

-  (2t  -  /So  «o)  (/Scoo)  cos  (/So  «o)l  t/  (2T/tDo  -  /So  -  |8) 

+  (i  x/a>o)I  (n  +  1)  J„  (/3<i>o)  sin  (/Sq  <oo)  -  /3  a»o  i„  _i  (/8  wq)  sin  (/So  wo) 

-  Pqo}oJ„  (/S  coq)  cos  (/So  <oo)1  U  (0q  -  0) 

+  (iT/2wo)l  -(n  +1)J„  (0too)  sin  Oo  coq)  +  /8  wo  ■^n-i  03  wo)  sin  (|8o  uq) 

+  03o  -  2t)  7„  (/S  «o)  cos  (00  <»>o)l  ^  03o  ~  2x/wo  -  3) 

+  (it/2o>o)  (-(«  +  T)  (0  «o)  sin  (/So  a>o)  +  0  '4>o  -^n  -1  (/3  (^)  sin  (0o 

+  (2x  +  00  «o)  -4  (/3  «o)  cos  (00  ««>o)l  V  (2x/ci>o  +  0o  -  0) 
when  n  =  0,  or  even  integer, 

=  (x/2a;o)  [-  (h  +1)  ^  (|3«o)  cos  (0o  wo)  +  0  «o  ^-i  (/3«o)  cos  (0o  tuo) 

+  (2t  -  00  «o)  Jn  03«o)  sin  (0o  «o)J  U  (Z-wIwq  -  0o  -  0) 

+  (x/wq)  [-(/i  +1)  J„  (0coo)  cos  (00  uo)  +  0  Wo  ■/«-!  (/3  wo)  cos  (00  Wo) 

—  00  Wo  Jn  (/3  wq)  sin  (0o  wo)]  U  (0©  —  0) 

+  (x /2ci>o)  [(n  +  1)  Jn  (0  Wo)  cos  (0o  Wq)  —  0  Uq  J„-i  (0  Wq)  cos  (00  Wq) 

~  (00  Wo  —  2x)  y„  (0  Wo)  sin  (0o  wo)J  U  (0o  —  Ivtuo  —  0) 

+  (x/2wo)[(n  + 1)  "/fi  (0  Wo)  cos  (0o  wq)  —  0  wq  J„-\  (0  wo)  cos  (0o  wo) 

+  (2x  +  /So  Wo)  (0  Wo)  sin  (0o  wo)]  U  (2x/wo  +  0o  ~  0) 


(55a) 


(55b) 
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when  n  =  odd  integer. 


ei,c  J  I  ^c(w)  \  ^  Jn  (fi  sin  (/So  ui)  d  W 

—  Qp 

=  (it/2wo)[— («  —  1) ./«  (^  <!))  sin  (/So  ci>o)  +  /3  wo  ^n-i  (/Swo)  sin  (/So  «o) 

-  (2»  -  /So  (do)  03  «o)  cos  (/So  (do)l  (2t/uo  -  /3o  -  /3) 

+  (iT/(do)  I(n  - 1)  ^„  (^  (do)  sin  (/So  (do)  -  /S  (do  -1  (0uo)  sin  (/So  (do) 

-  00  030  J„  (/S  (do)  cos  (/So  (do)]  V  (jSo  -  /3) 

+  (iir/2(do)[—  (w  ~1)  •/«  (^  (do)  sin  (^o  “o)  +  /3  (do  7„_i  (0o3q)  sin  (0q  (do) 

+  (00  (oo  -  2t)  J„  (/S  (do)  cos  (/So  (do)]  1/  (00  -  2T/(do  -  0) 

+  (ix/2(do)[-  (n  -l>/„  (0  (do)  sin  (/So  (do)  +  0  (do  J„.i  (jS  (do)  sin  (/So  (do) 

+  (2t  +  ^0  («>o)  Jn  (0  «o)  cos  (00  (do)]  t/  (2T/(do  +00“  0), 
when  n  =  0,  or  even  integer 

•^3.c  =  (*’/2(do)[  —  (n  —1)  y„  (0(do)  cos  (^0  t^)  +  00  Wo  "^n-I  (0  “o)  COS  (/So  (do) 
+  (2t  -  /So  (do)  J„  (0  (do)  sin  (/So  (do)]  U  (2x/(do  -  0o  -  0) 

+  (x/(do)[-(n  - 1)  y„  (0  (do)  cos  (/So  (dol  +  00  (do  /„  _1  (0(do)  COS  (/So  (do) 

—  00  (do  Jn  (0(do)  sin  (/So  (do)]  U  (/So  —  |3) 

+  (x/2(do)  I(«  -  l)y„  (0  (do)  cos  (/So  (do)  -  0(do  _1  (0  (do)  cos  (/So  (do) 

~  (00  ((Jo  ~  2x)  J„  OS  (do)  sin  (/So  (do)]  U  (/So  -  2x/(do  -  0) 

+  (x/2(do)[(/I  - 1)  J„  (/3(do)  cos  (/So  (do)  -  0  (do  y„  _1  (/S  (do)  cos  (00  (do) 


(56a) 
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+  (2t  +  /S  wo)  03  a>o)  sin  (/Sq  wq)]  U  {2t/wo  +  ~  P) 


(56b) 


when  n  =  odd  integer. 
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a  =  z/(d2/Xq) 


Fig.  3a  —  Normalized  energy  densities  along  the  axis  of  the  disk  for  a  <  1 


CW  FAR  FIELD  BEGINS 


a  =  z/(D^/Xo) 

Fig.  3b  —  Normalized  energy  densities  along  the  axis  of  the  disk  for  a  >  1 


ORMALIZED  ENERGY  DENSITY 
FOR  SINE  PULSE 


a  =  z/(D2/Xq) 


Fig.  4a  —  Normalized  energy  densities  associated  with  the  single-cycle  sinwo^  fo*"  different  values  of  s 


NORMALIZED  ENERGY  DENSITY 
FOR  COSINE  PULSE 


a  =  z/(D^/>,q) 


Fig.  4b  —  Normalized  energy  densities  associated  with  the  single-cycle  cos«o^  for  different  values  of  s 
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NORMALIZED  FIELD  FOR  SINE  PULSE 


Fig.  Sa  —  Normalized  field  as  a  function  of  r  for  a  =  0.2S 
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NORMALIZED  FIELD  FOR  SINE  PULSE 
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T  =  in’  =  NORMALIZED  TIME 

Fig.  5b  —  Normalized  field  as  a  function  of  t  for  a  =  0.5 
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NORMALIZED  FIELD  FOR  SINE  PULSE 


0=1 


Fig.  Sc  —  Normalized  field  as  a  function  of  r  for  a  =  I.O 


52 


